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The bouncing universe provides a possible solution to the Big Bang singularity problem. In this 
paper we study the bouncing solution in the universe dominated by the Quintom matter with an 
equation of state (EoS) crossing the cosmological constant boundary. We will show explicitly the 
analytical and numerical bouncing solutions in three types of models for the Quintom matter with 
an phenomenological EoS, the two scalar fields and a scalar field with a modified Born-Infeld action. 

S' 

O I. INTRODUCTION 

o : 

A bouncing universe with an initial contraction to a non- vanishing minimal radius, then subsequent an expanding 
phase provides a possible solution to the singularity problem of the standard Big Bang cosmology. For a successful 
bounce, it can be shown that within the framework of the standard 4-dimensional Friedmann-Robertson- Walker 
(FRW) cosmology with Einstein gravity the null energy condition (NEC) is violated for a period of time around the 
bouncing point. Moreover, for the universe entering into the hot Big Bang era after the bouncing, the EoS of the 
matter content w in the universe must transit from w < — 1 to w > - 1. 
' ' ' The Quintom model [l[ , proposed to understand the behavior of dark energy with an EoS of w > — 1 in the past 
O -1 and w < — 1 at present, has been supported by the observational dataQ. Quintom is a dynamical model of dark 
energy. It differs from the cosmological constant, Quintessence, Phantom, K-essence and so on in the determination 
of the cosmological evolution. A salient feature of the Quintom model is that its EoS can smoothly cross over w = — 1. 
I— ~~ '■ In the recent years there has been a lot of proposals for the Quintom-like models in the literature. In this paper we 
I | study the bouncing solution in the universe dominated by the Quintom matter and working with three specific models 
. we will show explicitly the analytical and numerical solutions of the bounce. 

We will start with a detailed examination on the necessary conditions required for a successful bounce. During the 
ON ■ contracting phase, the scale factor a(t) is decreasing, i.e., a(t) < 0, and in the expanding phase we have a(t) > 0. 
At the bouncing point, a(t) = 0, and around this point d(t) > for a period of time. Equivalently in the bouncing 
cosmology the hubble parameter H runs across zero from H < to H > and H = at the bouncing point. A 
■ successful bounce requires around this point. 

o : 

. H = -AnGp(l + w) > . (1) 

o ■ 

I From ([T]) one can see that w < — 1 in a neighborhood of the bouncing point. 
• , After the bounce the universe needs to enter into the hot Big Bang era, otherwise the universe filled with the matter 
with an EoS w < — 1 will reach the big rip singularity as what happens to the Phantom dark energy Q. This requires 
^ , the EoS of the matter to transit from w < — 1 to w > - 1. 

In this paper, we study the bouncing solutions in the Quintom models. The paper is organized as follows. In section 
II, we present the analytical and numerical solutions for different types of models of the Quintom matter. Specifically 
we consider three models: i)a phenomenological Quintom fluid with a parameterized EoS crossing the cosmological 
constant boundary; ii)the two-field models of Quintom matter with one being the quintessence-like scalar and another 
the phantom-like scalar; iii) a single scalar with a Born-Infeld type action. Ill is the summary of the paper. 
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II. BOUNCING SOLUTION IN THE PRESENCE OF QUINTOM MATTER 
A. A phenomenological Quintom model 

We start with a study on the possibility of obtaining the bouncing solution in a phenomenological Quintom matter 
described by the following EoS: 

w(t) = -r - - . (2) 

In @ r and s are parameters and we require that r < 1 and s > 0. One can see from ^ that w runs from negative 

infinity at t = to the cosmological constant boundary at t = \J arL d then crosses this boundary. 

Assuming that the universe is dominated by the matter with the EoS given by ([2|) , we solve the Friedmann equation 
and obtain the corresponding evolution of hubble parameter H(t) and scale factor a(t) as follows, 

H ® = In hi— ' ( 3 ) 

3(1— r)i + s 

a(t) = (f + -^)^J . (4) 

1 — r 

Here we choose t — as the bouncing point and normalize a = 1 at this point. One can see that our solution provides 
a picture of the universe evolution with contracting for t < 0, and then bouncing at t = to the expanding phase for 
t > 0. In Fig. [T]we plot the evolution of the EoS, the hubble parameter and the scale factor. 




FIG. 1: Plot of the evolution of the EoS w, the hubble parameter H and the scale factor a as a function of the cosmic time t. 
Here in the numerical calculation we have taken r = 0.6 and s — 1. 



One can see from Fig. [T] that a non-singular bouncing happens at t = with the hubble parameter H running 
across zero and a minimal non-vanishing scale factor a. At the bouncing point w approaches negative infinity. 



3 



B. Two-field Quintom model 



Having presented the bouncing solution with the phenomenological Quintom matter, we now study the bounce in 
the scalar field models of Quintom matter. However it is not easy to build a Quintom model theoretically. The No-Go 
theorem proven in Ref. [3 (also see Ref. 0,1, 1,0, HQ) 

forbids the traditional scalar field model with a lagrangian 
of general form C = C(<f>, V^V^) to have its EoS cross over the cosmological constant boundary. Therefore, to 
realize a viable Quintom field model in the framework of Einstein's gravity theory, it needs to introduce extra degree 
of freedom to the conventional theory with a single scalar field. The simplest Quintom model involves two scalars with 
one being the Quintessence- like and another the Phantom- like [l|, [l(| • This model has been studied in detail later 
on in the literature. In the recent years there have been a lot of activities in the theoretical study on Quintom-likc 
models such as a single scalar with high-derivative [HI, E2], vector field [HI, extended theory of gravity [lj| and so on, 
see e.g. [l5j . 

In this section we consider a two-field Quintom model with the action given by 



S = 



2<9 M <M M 2 - V (fa, fc) 



(5) 



where the metric is in form of (+, — , — , — ). Here the field <pi has a canonical kinetic term, but <f>2 is a ghost field. In 
the framework of FRW cosmology, we can easily obtain the energy density and the pressure of this model, 



and the Einstein equations are given by 



- -44 - V 



H2 = SjrGl 
3 V 



3^i + ^- = , 
d(pi 

dV 

- 3#0 2 - — = 
d(f>2 



(6) 

(7) 
(8) 
(9) 



From Eq. |T]), we can see that a bouncing solution requires </>?, = (j>\ + 2V when H crosses zero; and the Quintom 
behavior requires 4>\ = <j>f when w crosses — 1. These constraints can be easily satisfied in the parameter space of this 
model. 

In Fig. [2] and Fig. [3l we show the bouncing solution for two different type of potentials. In Fig. [2] we take 



= V x e 
<pi and 



-Ai 



+ V 2 e 



-a 2 



In the numerical calculation we normalize the dimensional parameters such 



62 by a mass scale M which we take specifically to be 10 M p i. And the hubble parameter is 



as Vi, V 2 , 

normalized with j^—- One can see from this figure the non-singular behavior of the Hubble parameter and the scale 
factor for a bounce. The EoS w crosses over the w = — 1 and approaches to negative infinity at the bouncing point. 
And due to the oscillatory behavior of the field <\>\ in the evolution, the EoS w is also oscillating around bouncing 
point. 

In Fig. [31 we take V{4>) — |m 2 0i 2 + Vo</>2~ 2 - This model also provides a bouncing solution, however the detailed 
evolution of the universe differs from the one shown in Fig. [2] Fig. [3] shows that the EoS of the Quintom matter will 
approach w = 1 asymptotically. 



C. A single scalar with high-derivative terms 



In this section we consider a class of Quintom models described by an effective lagrangian with higher derivative 
operators. Starting with a canonical scalar field with the lagrangian £ = ^d^ifid^cf) — V(<p). This type of models 
has been considered as a candidate for dark energy, however as shown by the No-Go theorem it does not give w 
crossing -1. As an effective theory we know that the lagrangian should include more operators, especially if these 
operators involve the term □</>, as pointed in Ref. [ll[ it will give rise to an EoS across w = —1. A connection of 
this type of Quintom theory to the string theory has been considered in Ref. [l6| and [T3] ■ In this paper we take the 
string-inspired model in [la ] for the detailed study on the bouncing solution, where the action is given by 



(10) 



s = 



XV -9 



-V(4>) vO^V^VWp 
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t t t 



FIG. 2: The plots of the evolutions of the EoS w, hubble parameter H and the scale factor a. In the numerical calculation 

we choose V((/>i,4>2) = Vie + y 2 e with parameters: Vi — 15, Vi = 1, Ai = —1.0, A2 = 1.0, and for the initial 

conditions <^>i = 0.5, cj>i = 0.1, (j>2 = 0.3, <^2 = 4. 




t t t 



FIG. 3: The same plots as Fig. [2]with different potential and model parameters V(cj>) = \rr?(\>\ 2 + V (f>2- 2 , m = 2, V = 0.4, 
and for the initial conditions 4>i =2, <j>i = 3, <f>% = 1, <f>2 = 2. 

This is a generalized version of "Born-Infeld" actionp"?], with the introduction of the (3' term. To the lowest 
order, the Box-operator term <f>\2<j) is equivalent to the term V M (/)V M when the tachyon is on the top of its potential. 
However when the tachyon rolls down from the top of the potential, these two terms exhibit different dynamical 
behavior. The two parameters a 1 and /?' in (|10p could be arbitrary in the case of the background flux being turned 
on [l9| . One interesting feature of this model is that it provides the possibility of its EoS w running across the 
cosmological constant boundary. In the analytical and numerical studies below to make two parameters (a 1 , j3') 
dimensionless, it is convenient to redefine a — a'M 4 and (3 — (3'M 4 where M is an energy scale of the effective theory 
of tachyon. 
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From |[T0|) we obtain the equation of motion for the scalar field 



^□(_r ) + aVAI (— p^) + M%/ + ^rD0 = O, (11) 



where / = ^/l — a'V AI 0V^</> + P'(pO<p and = dV/d<j). Correspondingly, the stress energy tensor of the model is 
given by 

T» v = g^[Vf - ^4 Vp(^V^)] + ^ + ^A_V M (^)V^ + V„(^)V^ . (12) 

Technically, it is very useful to define a parameter = = — ^u^f *° sorve CD ant ^ C3- In the framework of a 
flat FRW universe filled with a homogenous scalar field <fi, we have the equations of motion in form of 



Moreover, the energy density and the pressure of this field can be written as 



'=-/T -^-lm v - (15) 

P= ~^ ~ #+ 4M^ K (16) 



From Eq.fl]), one can see that a successful bounce requires: 



4M 4 V> 



V* + ^ = _^_^<0. (17) 



We will show below that (| 1T[) can be satisfied easily for our model. In the numerical study on the bouncing solution, 
we constrain the parameters a and (3 so that the model in (|10[) when expanding the derivative terms in the square 
root to the lowest order gives rise to a canonical kinetic term for the scalar field <fi [l6T |. i.e., a + (3 > 0. 

In Fig. 0] and Fig. [5] we show the bounce solution for different potentials. In Fig. [5J we take V(<fi) = Voe~ x ^ l M 
with A being a dimensionless parameter. One can see from this figure the scale factor initially decreases, then passes 
through its minimum and increases. Moreover, away from the bouncing point in the expanding phase, the EoS of 
the scalar field crosses w = — 1 and approaches w — —0.6, which gives rise to a possible inflationary phase after the 
bouncing. In the numerical calculation we take the energy scale M to be 10 -2 M p i, and the hubble parameter is 
normalized with M-. 

M p i 

In Fig. El we consider the model with potential V(<p) — Vo/(f) and then show another example of the bouncing 
solution 1 . Here the energy scale M is chosen to be \Q~ 2 M p i as well. One can see from this figure the clear picture of 
the bouncing, however the detailed evolution of the universe differs from the one shown in Fig. |4j After entering the 
expanding phase, the EoS w crosses the cosmological constant boundary and approaches w = |, which is equivalent 
to the EoS of the radiation. 



III. CONCLUSION AND DISCUSSIONS 



In this paper we have studied the possibility of obtaining a non-singular bounce in the presence of the Quintom 
matter. In the literature there have been a lot of efforts in construc ting the bouncing universe, for instance, the Pre 
Big Bang scenario (2l|. and the Ekpyrotic scenario [22j. In Refs. [23l \24. l25j and [2^, [23] the authors have considered 
models with the modifications of gravity with the high order terms. In general these models modify the 4-dimensional 



1 A Born-Infeld lagrangian with this potential provides a scaling solution, see Ref. [2dH . 
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FIG. 4: The plots of the evolution of the EoS w, the hubble parameter H and the scale factor a. Here in the numerical 
calculation we take the potential V(<p) = Vbe _A * > , a = —0.2, (3 = 2, A = 2, Vo = 5, and the initial values are <f> — 1, 4> — 3, 
H = -1, and ip = -80. 




FIG. 5: The plots of the evolutions of the EoS w, hubble parameter H and the scale factor a. In the numerical calculation we 
choose the potential as V(4>) = a = —0.2, (3 = 2, Vq = 0.7, and for the initial conditions cj> = 10, = —3, H = —1, ^ = —40. 



Einstein gravity. However, the models we consider for bounce universe in this paper are restricted to be within the 
standard 4-dimensional FRW framework. 

Recently two papers [111 [29| have studied the possibilities of having a bounce universe with the ghost condensate. 
In the original formulation the ghost condensate [30( will not be able to give EoS crossing w = — 1. The authors of 
these papers [1^, [2(| have considered a generalized model of ghost condensate [3l| and shown the bouncing solutions. 
In this paper we have studied the general issue of obtaining a bouncing universe with the Quintom matter. Our 
results show that a universe in the presence of the Quintom matter will avoid the problem of the Big Bang singularity. 
Explicitly for the analytical and numerical studies we have considered three models: the phenomenological model, 
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the two-field model and the string-inspired Quintom model. The latter one is a generalization of the idea in Ref.[ll| 
by introducing higher derivative terms to realize the EoS crossing w — — 1. In this regard, this model for the bounce 
solution has the similarity with a recent paper 32] where the authors presented a bouncing solution with non-local 
SFTfU. 
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